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We study multi-field inflation models that contain a non-trivial field-space metric and a non- 
minimal coupling between the gravity and infiaton sectors. In such models it is known that even 
in the absence of explicit interaction terms the infiaton sector can decay into matter as a result 
of its non-minimal coupling to gravity, thereby reheating the Universe gravitationally. Using the 
Bogoliubov approach we evaluate the gravitational decay rates of the infiaton fields into both scalars 
and fermions, and analyse the reheating dynamics. We also discuss how the interpretation of 
the reheating dynamics differs in the so-called Jordan and Einstein frames, highlighting that the 
calculation of the Bogoliubov coefficients is independent of the frame in which one starts. 


I. INTRODUCTION 

An epoch of inflation in the very early Universe is now firmly supported by recent observations of the cosmic 
microwave background (CMB) These observations suggest that the primordial curvature perturbation 

generated during inflation is nearly Gaussian and adiabatic, with a power spectrum that deviates from scale- 
invariance at the 5(T-level. They also suggest that the amplitude of tensor modes is relatively small, i.e. that 
the energy scale of inflation is low. 

With inflation widely accepted as a key part of the standard model of cosmology, the question now turns 
to determining the exact nature of inflation and how it might be embedded in some fundamental high- 
energy-physics theory. Inflation models containing non-minimal gravitational coupling comprise one class 
of promising models. As well as being theoretically well motivated in the context of high-energy-physics 
theories such as string theory, see e.g. 0, their predictions also lie at the sweet spot of current observational 
constraints Examples include Starobinsky’s original inflation (written in its scalar-tensor form) Q, 
Higgs inflation Q and a whole class of so-called conformal inflation models recently proposed by Kallosh et 
al. Q. Whilst most of these models are studied as single-field models, the high-energy-physics theories that 
motivate them generically predict the presence of multiple fields during inflation. As such, it is important to 
determine any possible signatures of multi-field effects in models with non-minimal coupling [l0l-[l5j. 

It has recently been demonstrated that in constraining specific models of inflation with current CMB 
data, details of the reheating process must be properly taken into account, even for single-field models 
|lfil - ll9l| . This is testimony to the precision of current CMB data. Moreover, in the context of multi-field 
models of inflation, the primordial curvature perturbation may continue to evolve during reheating, and this 
evolution must therefore be tracked until an adiabatic limit is reached [20l - l^ . It is known that reheating can 
take place gravitationally in models with non-minimal gravitational coupling; even if there are no explicit 
interaction terms between the infiaton sector and matter, gravitational particle production takes place as a 
result of the non-minimal gravitational coupling [3, 1234)^ . In light of the renewed interest in this class of 
models, in this paper we revisit the theory of gravitational reheating after inflation and present a multi-field 
formulation of gravitational particle production. Whilst we focus on perturbative gravitational reheating, 
we nevertheless employ the method of Bogoliubov transformations to determine the decay rates. As such, 
many of the features we discuss should also carry over to the case of non-perturbative preheating. Of course, 
the Bogoliubov approach recovers the standard perturbative quantum-field-theory (QFT) results when the 
appropriate limits are taken. 

Throughout the analysis we try to pay particular attention to how the interpretation of the reheating 
dynamics differs in the so-called Jordan and Einstein frames. During the oscillatory phase at the end of 
inflation, the Hubble rate in the original Jordan frame contains an oscillatory component, and it is this 
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oscillatory component that gives rise to particle production even in the absence of direct couplings between 
the inflaton sector and the decay products. The evolution of the scale factor in the Einstein frame, on 
the other hand, is equivalent to that of a matter-dominated universe, and can essentially be neglected. As 
such, the leading-order contribution to the gravitational particle production is not a result of the oscillatory 
nature of the Hubble rate. In its place, however, one obtains explicit gravitationally induced interaction 
terms between the inflaton sector and ordinary matter, through which reheating proceeds. Although the 
interpretation in the two frames is different, we nevertheless find that the calculation of the Bogoliubov 
coefficients is independent of the frame in which we start; working in conformal coordinates and requiring 
that the mode functions under consideration be canonically normalised leads us to a common set of variables 
and form of action. 

A technical complication that arises in the context of multi-field models with non-minimal coupling is that 
even if one starts with a flat field space in the Jordan frame - i.e. a canonical, diagonal kinetic term - then 
one obtains a non-flat field space in the Einstein frame, where evaluation of the inflaton dynamics is simpler 
[ 3 ^. We thus find it necessary to work in the mass eigen-basis as defined with respect to the Einstein frame 
potential. As a result, if either of the Jordan frame field-space metric or non-minimal coupling are functions 
of some light spectator field, we find that t he g ravitational decay rates generically become modulated, giving 
rise to a modulated-reheating scenario 

The rest of this paper is organised as follows: In Sec.|TT]we begin by outlining the class of models under 
consideration and by reviewing some of their key characteristics, with the review extending into AnnendixIAl 
In Sec. uni we then analyse the reheating dynamics. We start, in Sec. IIII A[ by looking at the background 
dynamics of the oscillating inflaton fields at the end of inflation, and in Sec. IIII Bl and Sec. IIII Cl we turn to 
the reheating process itself, presenting the details of the Bogoliubov calculation used to determine the decay 
rates. Additional details regarding the calculation of fermion production rates are included in AppendixlBl 
Finally, Sec. lIVI is devoted to summary and conclusions. 


II. MULTI-FIELD MODELS WITH NON-MINIMAL COUPLING 

In this section we define more explicitly the class of models under consideration and also discuss the relation 
between formulations made in the Jordan and Einstein frames. 


A. Actions in the Jordan and Einstein frames 

The general class of models that we are considering take an action of the form 

S = J - \habg^'^d^rd.4>^ - + . 5 ™, (I) 

where a,b = I...n label n scalar fields that are potentially all non-minimally coupled to the Ricci scalar R 
through the function /(</>). hab defines a non-flat field-space metric and V is some general potential depending 
on all the fields. We take the matter part of the action to consist of bosons and fermions, namely 

Sxi = I - U{xz)} , 

S'™ = where f-^ - I 

I t 5 ^, = -■ 

Here 0 is given as ij) = with = ,9^ -p T^ and where is the tetrad defining 

local Lorentzian coordinates, 7 “ are the standard Dirac matrices satisfying { 7 “, 7 ^} = 277 “^, and the spinor 
connection T^ is defined as T^ = (I/ 2 )E“^e^V^e/ 3 A, where = ^[ 7 “, 7 ^]. We also have = ipiP, where 
/? = We have omitted the conformally invariant gauge fields, as they do not play an important role in 
the perturbative reheating considered in this paper. See, however, [3l| for a discussion on the gauge trace 
anomaly and its importance in the reheating process. 


^ Note that we are working with the signature (—h H—h). 
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Matter in the above action is minimally coupled to gravity, and this “frame” is referred to as the Jordan 
frame. However, on making the conformal transformation g^i, = with = /((/>)/Mpj, the action can 

be re-written as 


S = J 


j -vl+SrH, 


2 2 


whercQ 


_ -^pi f, ^fafb 

^ab — r I ^ab \ 


2 / 7 ’ 


and the matter actions now take the form 


P 

u{x^y 


5"^, = - y 0 tpi -f ■ 


Here we have defined 

= = + Xidpinil) and 0 = 6 ^ 7 “ (3^-k T^), 


X^ 


n 


( 3 ) 


( 4 ) 

( 5 ) 

( 6 ) 

( 7 ) 


where which gives x^{x) = 7 '^(x)/H, and the spinor connection is conformally invariant (see, 

e.g., footnote 4 of [Slj). In this form, the fields (/)“ are minimally coupled to gravity and the gravity sector 
is of the standard Einstein-Hilbert form. The matter sector, however, becomes explicitly coupled to the 
inflaton sector, and we must also be careful to take into account the spacetime-dependent rescaling of units 
that results from the conformal rescaling of the metric. 


B. Einstein’s equations and the equations of motion 


Having defined our actions, let us briefly review the gravitational equations of motion that they give rise 
to. In this section we will simply quote the main results. Additional details regarding these known results 
can be found in Appendix [Al 
Re-expressing (HD in the form 




where contains the matter sector, and minimising ([SI) with respect to g^'^ we get 

I 


G,u - ^ L 


Tiyj+P:;^+X7^V.f-g^.af 


where 


7^(b) = 

in' 


2 




Similarly, varying the action with respect to the fields P we get the equations of motion 

habOp + r,„l,gf^’'VppV,P -Va+ faR = 0, 


( 8 ) 


(9) 


( 10 ) 


( 11 ) 


where Tab\c = Tfec i® Ih® Christoffel connection associated with the field-space metric hab- 


Note that, for example, fa denotes taking the derivative of / with respect to the a’th field. 


2 
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Turning to the Einstein frame, we can similarly write the action ([3]) in the form 

5 = + 5^^ =-^Sabr''d^,rd.(IJ^-V, ( 12 ) 

and we will return shortly to the relation between and We then find the standard Einstein 

equations 



where 


'fW 

fiiy 



Sab^u(l)^ - Otiu 


Sabr’^Vpr^o 



(14) 


The equations of motion for the fields in the Einstein frame take the form 

- SabUctb^ + 14 + W = 0, (15) 

where ^'®^rf,c|o = and is the Christoffel connection associated with Sab- 

Regarding the matter energy-momentum tensors, one can show (see Appendix A for a review) that under 
certain conditions the following relations hold: 

= V'"tM=0 and (16) 

For future reference we note that the energy-momentum tensors associated with Xi s-iid ipi are given, 
respectively, as 


+ U{xi)^ , (17) 

^ (V'i7(/.(2;)-D,.)V'i - , (18) 

where our symmetrisation with respect to the indices includes a factor of a half. In the Einstein frame we 
similarly have 


= 'DpX^V,X^ - 9p. {^r^VpX^VaX^ + , (19) 

^ {i’iX{p{x)D^)'^i - D(^^'4}i%){x)'ipi^ . ( 20 ) 

III. REHEATING DYNAMICS 

Having described the general setup for our class of models and the relation between the Jordan and Einstein 
frame formulations, in this section we turn to the process of reheating. We will begin by considering the 
background dynamics of the inflaton fields after the end of inflation, before then turning to the particle 
production process. We also consider the effect of the produced particles on the dynamics of the inflaton 
helds and how reheating ends. At every step we try to discuss how the interpretation of the reheating 
process differs in the Jordan and Einstein frames. For a review of reheating after inflation and the techniques 
employed in this section, see e.g. [l3, IdOl. 1^. 


A. The oscillating phase 


In describing the dynamics of the inflaton fields after the end of inflation we make the standard assumption 
that at background level our Universe is described by a Friedmann-Lemaitre-Robertson-Walker (FLRW) 
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metric. Indeed, the homogeneity and isotropy of the Universe should be guaranteed thanks to the preceding 
epoch of inflation. If we would like to write both the Jordan and Einstein frame metrics in FLRW form, then 
the relation g^i, = gives us 

ds^ = —dP + cP{i)5ijdx^dff = = VP {—dP + P{t)5ijdx'‘dx^) , (21) 

where we recall that = /(</))/Mpj. From the above equation we then find the following relations: 

a = Uo, dt = Qdt, dP = dx^ and ^ ~ ^ 

If we consider the epoch before reheating, when only the inflaton sector is present, then the Friedmann 
equation and equations of motion for the scalar fields in the Jordan frame are given, respectively, as 


= 


+ V-3Hf 


and 


dt 


+ (Vb - fbR) = 0 , 


(23) 


where /dt = (j>°‘ + and a dot denotes a derivative with respect to the Jordan frame cosmic time 

t. Similarly, in the Einstein frame we have 


3iJ^ = 


Ml, 


1 dcf)'^ d(p 
2^°'''~£~dd 


■V 


and 


D{d<j)°^/di) 
di 


3H 


d(j)°‘ 

di 


+ S^'^Vb = 0 , 


(24) 


where D{d(j)°-/di)/di = d'^(j)°'/dP + /di){d4>‘^/di), V = M|,U//^ and i denotes the cosmic time in 

the Einstein frame. 

Given that the non-minimal coupling between the inflaton fields and gravity is removed in transforming 
to the Einstein frame, it is much more convenient to solve for the inflaton dynamics in this frame. As such, 
let us proceed by first solving for the dynamics in the Einstein frame. Do note, however, that it should be 
possible to solve directly in the Jordan frame, see e.g. [s^. 

Given that in the Einstein frame the dynamics are determined by the Einstein frame potential V and 
the field-space curvature associated with Sab, "we make the assumption that at the end of inflation all of 
the inflaton fields begin to oscillate about the minimum of V at cj)°- = and can be decomposed as 

(j)°- = -\- cr“. Requiring the absence of a cosmological constant dictates that Mev = 0. Gombining this 

with the fact that Va\vev = 0, on expanding the inflaton part of the Einstein-frame action to second order in 
cr“ we get 


5* = 



M|iR 

2 


7:Sab\vev9^''dfj,a°-d^a^ - it4b|vevCT“cr'’ 


(25) 


where we have made the assumption that the potential can be well approximated as being quadratic about its 
minimum. In order to deal with the non-diagonal nature of this action, we now introduce the mass eigenstates 
of the Einstein-frame potential. Namely, we take fT“ = , where 

= (26) 

with U^blvev = 5'“'^|vevl4b|vev and 5 ab|veve( 4 e^ = Jab- In the above expression the hat on the index A 
suppresses summation. The Einstein frame action then takes the form 

so that the equations of motion for oP are simply given as 



Note that in deriving the above results we have taken cr“ and hence to be background quantities that 
only depend on time, i.e. the decomposition </>“ = -|- (t“ is not a decomposition of into its classical 


da^ 

di 


A\2 


-mi(a-) 


(27) 
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background part and quantum perturbation, but is simply an expansion of the classical part of about (/>vev 
Also note that no terms involving the field-space curvature appear in (1281) . This is because such terms would 
be second order in (t“, taking the form /di)[d(j'^/£). As such, we see that models with very 

large field-space curvature, such as that considered in Fig. 5 of [1^, lie beyond the scope of our perturbative 
approach. 

Making the assumption m\ ^ dH/dt, i.e. that the timescale of the field oscillations is much shorter 

than that of the background evolution of the universe, we find the solutions 


cos[m^t -I- d^]. 


(29) 


where dA are constant phases. The Einstein frame Friedmann equation then gives us 


= 




E 


da^ 

di 


A\2 






A 


6Ml 


1 + 


pr 


3id 

2mA 


sin(2(TO^t -I- dA)) + 0{H^/w?a) 


(30) 

As such, we see that to leading order in H/mA the evolution of the Einstein frame Hubble rate coincides 
with that of a matter-dominated universe. On calculating dH/dt one finds 


did 

di 



1 


E 

A 




cos(2(m.4t -I- dA)) 


'pr 


(31) 


The second term in the brackets represents an 0(1) deviation from the case of matter-domination, which can 
be seen by noting from (13(1 that (oq )^m3^/(6Mpjd^) ~ 0{H^). However, the important result as far as we 
are concerned is that dH/dt ~ <C m\. 

We are now interested in using these results to determine the background evolution in the Jordan frame. 
As is evident from (Ea), in order to do this we need expressions for / and its derivatives, and these can be 
obtained by expanding / about /vev In doing so, we make the assumption that by the end of reheating, 
when all helds have decayed, /vev = -^pp We therefore have 


/ = M|i 



/a^° 


1 /abcr“cr'’ 

2 + 


= Af, 


PI 





1 fABCt^Ct^ 

2 


(32) 


where /a = fa^A- Inserting this expansion into the last relation in (12211 and evaluating to leading order in 
and H/mA we get 


If we assume that /a/Mpi ^ 0(1), and recall from (15(1 that a^mA/{Mp\a?/'^) ^ 0{H), we see that the 
evolution of the Hubble rate in the Jordan frame has an oscillatory component that is not suppressed. 
Note that to leading order in H/mA the cosmic times as defined in the Jordan and Einstein frames are 
interchangeable. With this, we see that H picks up a term that is 0{mAH) (assuming Ja/Mpi ~ 0(1))- 
This is to be compared with the case in the Einstein frame, where dH/di ^ 0{H^). 


B. Perturbative QFT approach to reheating 

Having discussed the background dynamics of the oscillating inflaton fields at the end of inflation, there 
are essentially two ways in which we can now consider reheating into ordinary matter. The first follows the 
standard perturbative QFT approach, and appears natural in the Einstein frame. The second method involves 
calculating Bogoliubov coefficients in an approach based on QFT in a time-varying classical background. This 
second method appears natural in whichever frame we begin, but the interpretation in each frame is somewhat 
different. In the case of perturbative reheating both methods are equally valid, and the result is independent 
of the method used. 
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Decay rates 


In transforming to the Einstein frame, one consequence of the conformal transformation is that we explicitly 
see the appearance of interaction terms between the inflaton sector and ordinary matter. These are apparent 
in the factors of that appear in 2?^, U/il'^ and in and Using the expansion of / given in 

dSll), and taking 


u{x) 


2 2 

2 


(34) 


we find that the Einstein frame action contains the tri-linear interaction terms 


A 


^int 


fAa 


(2 


mZ + 


and 


rv^ _ 

^'i'n + 


fAa 


2 M 2 j 


A _ _ 


(35) 


where we have integrated by parts and used the equations of motion for in deriving the first of theseH Here 
we neglect to consider four-point interaction terms, as we know that such terms cannot allow for complete 
reheating (See, however, Sec. IV of [^.1 

Another key feature of the Einstein frame is that the scale factor is evolving slowly, i.e. dH/dt <C m^, 

which allows us to neglect the expansion of the Universe. As such, we can use flat-space QFT calculations 
to determine the transition amplitudes for -A xx a"^ -A that result from the interaction terms 
in (ESI). These amplitudes can in turn be used to calculate the decay rates per unit time and volume of the 
oscillating fields 


where 


>xx 




1_4^ 


1/2 


and 


9xA — 


fae\{m\ + 2ml) 


4A/2j 




O^pA'^A 

Stt 


1 - 


4m; 

9 

m I 

A 


3/2 


and 


fae\m^ 


(36) 


(37) 


^pi 


In this approach we interpret the oscillating inflaton fields as a condensate of zero-momentum particles that 
can decay into two scalars or a fermion-anti-fermion pair. Our reason for suggesting that this approach seems 
“natural” in the Einstein frame is that it is in this frame that the necessary interaction terms are explicit 
and that the background evolution of the scale factor can be neglected. 


Dynamics including decay products 


Once the rate of decay becomes significant, namely Ta^H (see for the definition of f a), we must 
take into account the effect that the decay products have on the oscillating inflaton dynamics. Remaining 
in the Einstein frame, we see from m that the dynamics of the inflaton fields is sourced by the trace of 
the matter energy-momentum tensor. In Sec. IIII Al we ignored this term, assuming that inflaton decay was 
initially negligible, but now we must properly include it. At the level of the action the effect of matter fields 
on the dynamics of is evident in the explicit interaction terms, such as those given in (1351) . As such, in the 
context of the perturbative QET approach it is necessary to calculate 1-loop corrections to the propagator of 
a^. Invoking the optical theorem, one finds that the effective equations of motion for take the form [4fll | 



where 








(39) 


^ Note that as we have used the background equations of motion for the effective interaction term is only valid in making 
tree-level calculations. If we wish to go beyond tree-level calculations, then we would have to use the derivative interaction 
term directly. 
















On inserting the zeroth-order solution for H, namely H = 2/3t, the last term in the square brackets vanishes. 
If we also assume mA ^ r^, then the solutions to the above equations take the form 


. 3/2 


exp 


2 


cos \jn^t + . 


(40) 


Comparing with we see that there is an additional exponential decay of the amplitude of the oscillations. 

Phenomenologically, the effect of inflaton decay is often modeled by including an additional frictional term 
in the equations of motion for as follows (4nl| : 


d'^a^ 

di"^ 


(3i/ + f^) 


da^ 

di 


= 0 . 


(41) 


Indeed, under the assumptions ruA ^ H,Ta, one can see that & does satisfy this equation. The advantage 
of using this phenomenological equation is that it can be recast in a form that is intuitive. On multiplying 
through by da^/di and averaging over many cycles, it can be re-written as 


—h ^HpA + — 0, 

dt 

where we have once again assumed mA ^ Ha and 

I fda^V 1 


PA = t; 


2 \ dt 




Summing over all A we have 


where 


dp, 


dt 


J- + 3Hpa + TaPa — 0 , 


= '^PA and fg = ^Pa- 


(42) 


(43) 


(44) 


(45) 


We thus see that the energy density of the oscillating fields decays as a result of the Hubble expansion and 
the decay into matter particles, which is intuitively what we expect. It is important to note, however, that 
this phenomenological approach relies on the nature of the interaction terms considered and the fact that 
the inflaton fields are oscillating in a quadratic potential. As such, the situation will be different in the more 
general case [i^ . 

Assuming that the decay products quickly thermalise and can be modeled as a relativistic fluid, we also 
have the equations 


dPx. 

di 

dptpi 

di 


+ 4ffPx, 


(46) 

-|- AHp^. 


(47) 


A 


Combining with (I42p . one can see that the total energy density is thus covariantly conserved, in agreement 
with (fOl) . 


C. Bogoliubov approach to reheating 

In the flat-space perturbative QFT approach discussed in the previous subsection, we considered the 
oscillating inflaton fields as a collection of massive zero-momentum particles decaying into matter fields. 
One of the limitations of this approach is that it can only be applied in the perturbative regime, where 
interaction terms are small. An alternative approach to calculating the decay rates is based on QFT in a 
time-varying classical background. Within this framework particle production is a collective phenomenon, 
and the interaction terms do not necessarily have to be small. In the case that they are small, the results of the 
previous subsection are recovered, but in the case that the interaction terms are large it is possible to obtain 
resonant particle production, or preheating - see e.g. (4^. In this paper we will focus on the perturbative 
regime and confirm agreement with the perturbative QFT results given in the previous subsection. Much of 
our discussion, however, will also be relevant in the preheating regime. 
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Jordan and Einstein frame interpretations 


In calculating particle production in a time-varying classical background one is interested in solving for 
the mode functions of the matter field under consideration, and it is important that these are the 

mode functions associated with canonically normalised fields. As a consequence, we find that the calcula¬ 
tion becomes independent of the frame in which one starts, leaving only a difference in interpretation. To 
demonstrate this let us consider the case of the bosonic field %. 

Specialising to the case of an FLRW metric, the Jordan frame action for y becomes 


Sx= dtdrxa'^- 


1 




(48) 


where we have assumed U{x) = w^X^/2- In order to bring this into canonical form we use conformal time - 
defined as adrj = dt ~ and also introduce the re-scaled field u = ay, giving 







(49) 


where a prime denotes differentiation with respect to conformal time. Working in Fourier space this gives 
rise to the equations of motion for the mode functions as 

a” 

u'l + wluk = 0 with wi=k‘^+ a^rn^ -, (50) 

^ a 

where k = \k\. As already discussed in Sec. lIII Al the scale factor in the Jordan frame has a rapidly oscillating 
component, and it is the resultant rapid time-variation of the effective mass of Uk that gives rise to particle 
production. 

If we now go to the Einstein frame, the action for x takes the form 


Sv= did^xd^ 


di 


1 


- - m^x^ + 


2M2j 


( 2 ^ 


m. 


2 \ ~2 
X 


(51) 


where we have expanded / to linear order in a^, integrated by parts and used the equations of motion for 
in order to get the interaction terms. In order to bring this into canonical form we once again use conformal 
time (note that conformal time is frame independent) and define the field u = a-x- However, given that 
X = x/H, we see that 


a 

u = ax = -^X = ax = u- (52) 

As such, whichever frame we start in, the analysis becomes identical once we transform to the canonically 
normalised variables. Indeed, on using (I22|) one finds 


= a 


(^H + 


2 i\/r2 
PI 




/ 

di 


Mpj 


H^\ - 


df 


1 (dy ^jjdf 

2 A/|jVdi 2 ' AfMl^\di 


(53) 


Expanding / to first order in and using the equations of motion for this reduces to 


a" d" fAa ~2 2 


a a 2Ml 


PI 




(54) 


so that on substituting into dSO]) we have 


2 7,2 I ~2 2 d fAOi 

Wj, ~ k + a - 1 ---pr 

^ a Mi 


1 


PI 


~2 2 ^ J ~2 2 


(55) 


and this is in agreement with the frequency we would obtain from (1511) on defining u = ay. 

The only difference between the two frames, therefore, is the interpretation. In the Jordan frame the 
effective mass of the scalar field is oscillating as a result of the oscillating scale factor, which is why we refer 
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to the process as gravitational reheating. In the Einstein frame, however, the scale factor is slowly varying 
and the oscillatory nature of the effective mass of the scalar simply results from the explicit interaction terms 
in the action. 

Whilst above we have considered the case of a scalar field, one can easily see that the same applies for 
fermions. On using conformal time, in the Jordan frame the canonically normalised field is 'h = In 

the Einstein frame, on the other hand, we have 





= = V]/. 


(56) 


The production of bosons 


We now turn to calculating the particle production rate for bosons, and we follow very closely the analyses 
given in [10, . Let us start by expanding the quantum operator u = ax in the standard way as 


u = 


(27r)3/2 


Uka^t 


^ik-: 


Ukaj-e 


(57) 


where Uk satisfy (1501) and at and are the creation and annihilation operators satisfying the standard 
commutation relations. 

The Hamiltonian associated with the action (091) can then be expanded as 


H, = 


1 

2 



+ 


+ Ffeara_r + F>la 




k —k 


(58) 


where 


2Ek = KP + wl\uk\^ and Fk = ( 4 )^ + 


(59) 


It is possible to construct mode functions that satisfy Fk = 0 - thus diagonalising the Hamiltonian - as 


/ ^ o.k{v) 

Uk{r]) = , exp 


\/ 2 wfeM 

where ak{r]) and /3kiv) must satisfy the equations 


drj'wkiv') 


Pkijl) 

\/2w/c(r?) 


exp 


dT]'wk{j]') 


(60) 


! r \ 

o^kijl) = ^T-^exp 
2wk 


Wi. 


Id'kiv) = exp 
2wk 


2i 

-2i 


— OO 

rv 


drj'wkir]' 

dij'wkiv') 


Pkiv), 

ak{ri). 


(61) 

(62) 


We also require that |afep— = 1 in order that the canonical commutation relations for u are satisfied. One 
can confirm that these mode functions also satisfy the equations of motion, and that Fk = Wk{i/2+\Pk{v)\'^)- 
If in the asymptotic past Wkirjo) is approximately constant, then the mode functions defined above are a 
linear combination of the positive- and negative-frequency mode functions associated with the Bunch-Davies 
vacuum. If Pkivo) = 0 at this time, then the mode functions do indeed coincide with those of the Bunch-Davies 
vacuum, where Fk is minimised. As time evolves, however, the evolution of Wk{r]) causes Pkiv) to evolve 
away from zero, meaning that Fk is no longer minimised. Given the diagonal nature of the Hamiltonian, this 
effect can be interpreted as particle production. 

Let us now assume that at the initial time rjo at which inflaton oscillations commenced there are no x 
particles present, i.e. Pkivo) = 0 ^-ud a.kir]o) = 1- We then assume that at times shortly after rjo we are 
in the perturbative regime where Pkifl) ^ 1 ^-ud OLk{rf) — 1 <C 1. This ensures that we are in the same 
perturbative regime for which the QET calculations of Sec. IHI HI are applicable, where Bose condensate 
effects are neglected. Solving the above relations iteratively, we find a solution for /3fc(r?) as 



-2i f dr]"wkir]'') 

J —GO 


( 63 ) 























(64) 

(65) 


Using (EH), we find that w'^ and wf. to leading order in and H jmA are given as 


yj' ~ fAa^' 

^ 2wk 2M|j 

2 7 2 , -2 2 
Wj^ — rC CL ■ 


{m\ + 2ml) + 


aP'Hml 

Wk 


Taking expressions to leading order in ensures that we are only considering the tri-linear interaction terms 
and perturbative regime appropriate for comparison with the perturbative QFT calculations of Sec. IIII Bl In 
making order-of-magnitude estimates, we note that from (1301) one can deduce the order-of-magnitude relations 
a'/Mpi ^ aH and ~ H. We have also assumed that fAlM-px ~ 0(1) and k ~ 0{amA)- The 

first of these assumptions is compatible with expansion (15^ so long as we have jM^x 1 , and the second 
assumption comes from our expectation that modes with k ^ 0{amA) will be produced. On substituting 
dH, (ESI) and (EH) into (IH!?1) we notice that in general the integrand is highly oscillatory in 77 . As the second 
term in (|64l) is non-oscillatory, we find that its contribution to averages to zero. The first term in (1641) . 

however, is oscillatory, which allows for the possibility of stationary points in the total phase of the integrand. 
We thus have 


o , ^ sr^mAfAa^i^ml + ml) p , ^ ^ 

= E- k^ + &^irj')ml )] " ^xp ^mA^Pk^^ )] } , (66) 


where 


f drj" 

J —00 

= - r dr," + 

J -00 V mAjm^ 


2wk\ , dA 

;-I H-, 

mA J 


(67) 

( 68 ) 


The phase tp'k 2 iv) does not have a stationary point for physical values of a, so only the term involving ipkiW) 
contributes to I3k(j,)- Using the stationary phase approximation we find 


A 

where /3^ = 


.A ( 2^2 +m|) 4 


27r 


16M|jf 




(69) 

(70) 


r,^ is the time at which dip-^^/dr, = 0 for some given k and is the sign of fpkiiVk)- phase is stationary 
when 2wk = dm a, which to leading order gives 


d?{r,] 


1_4^ 


(71) 


Given that k/d{r]^) coincides with the momentum of the produced particle, this result coincides with our 
expectation from kinematic considerations. The second derivative of the phase is given as 




ml ^ 2M|j ^ " xy 


(72) 


The two step functions in the above expression for jdk {r,) simply reflect the fact that a certain mode will only 
have been excited ii r, > r,^ > rjo. Note that r,^ is different for different A. 

In looking to determine the production rate of x particles let us start by considering the continuity 
equation for the Einstein frame energy-momentum tensor associated with y. Under the assumption of a 
FLRW background, we know that the expectation value of the energy-momentum tensor can be expressed 
as ( 0 |f«)M^| 0 ) = diag( —Pxy Px^ Pxy Px)’ where |0) is the vacuum state as defined with respect to g^:. The 
continuity equation m can then be written as 
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In the absence of particle production, i.e. in the absence of the interaction terms given in (I35L the right-hand 
side of this equation would be vanishing. As such, in order to determine the particle production rate, we 
wish to evaluate the right-hand side of (1751) . 

Using (O to determine re-expressing the result in terms of the canonically normalised field u and 

taking the vacuum expectation value, we find 


h (5+1*1") - 


-Px + 3px = - 


(27r)^ 


+ + 27^2) Q + |^^|2 ^ 


~(w' + ">U) (5 




,-2i / Wkdr)' 


(74) 


(75) 


where Ti. = a'/a, 5ft(A) denotes the real part of X and 3(^) similarly the imaginary part. Assuming that 
Q ^ is of the same order of magnitude as 5ft ^ , keeping terms only linear in 

Pk and neglecting the vacuum density contribution, we find that to lowest order in'H/wk the right-hand side 
of (ITHll is given as 


1 , o~ ^ 

2fdi^ Px + 3Px)- 


1 df 


(27r)3a4 2/ di 


j d^k^ [2wl+mld^] 5ft . 


(76) 


In general this quantity is highly oscillatory, and we are therefore interested in finding its average over several 
oscillations, i.e. over a time-scale T ~ C>(l/(aTO^)), where we have chosen to work in conformal time. Taking 
CTfc ~ 1 and substituting the results (17^ and (IMl) we have 


1 df , 
2fdP 


^Px 


1 r^+T' ^ 


47r 


2T 


V-T 


^ (27r)3d4 


dkk^Wk 


[2wl + mld^] fAajmA 

8Mpja3/2i 


(77) 


X 2 fS 


13: { 


^imAipk,lif) _ g'irnA'<IJk,2if) 


where () denotes taking the average over several oscillations. However, we can see that, due to the oscillatory 
nature of the integrand, this average will only give a non-zero result if r] coincides with a stationary point 
of the phase of one of the terms in the integrand. As there are no stationary points of '0^2 (ft) physical 
values of a, the only terms giving a non-zero contribution are those containing the phase 0^i(?7), namely we 
find 


2fdi 


{-Px + 3Px) ) = 


rnt+T 


47r 




rk Wk 


[2wl+mld^] fAa^rriA 


8M2jh3/2i 


(78) 


X 2 fS !(’?') 


Seeing as the integral over rj' is centred on the stationary point for each A, we can take T —>■ 00 , as contribu¬ 
tions away from the stationary point will average to zero. On making the stationary phase approximation, 
and after a little manipulation we eventually find 

(j^§(-ftx + 3Px)) =1] (27r)3l5(^) J dk6{7j-r]^)k^Wk{Vk)l3kl3k* (79) 

X 2cos {mAipk,i{Vk) - mBtl^kAiVk) + (4 “ sf )V4) ©(ftfc - ftf )0(ftf - fto), 


where we have used (15^ and (1701) . If we arrange that niB > tua for B > A, meaning that 0(?7^ — t]^) = I 
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only ioi B > A, then dzi) can be written as 


V di 


{-Px + ^Px)} =Y^ 


47r 


(27r)3a®(?7) 


E 

A,B>A 


47r 


(27r)3a^(?7) 


dk5{ri - r]^)k^Wk{rit)\l3^\'^ 

^ dkS{ri-7ji)k^Wk{Vk)PtPk* 


(80) 


X 2 cos {mAtpt,iiVk) - mBtpk,lidk) + (4 - sf )7r/4) ©(^if - Vo)- 

The delta function in r] can then be expressed as a delta function in k by using the fact that S{ri — r]^) = 
\'4’kiiPk)\^i'^k'iip)) combination with relation (ITTl) . and we hnd 


T~f~} ^ ^ 7T7 A i 

^{P-Pk) = Y^\'^k^{Vk)\^{k-d{n)^i^), where ^ 1 -^ 


2 \ 1/2 


(81) 


If we consider only the diagonal contributions to the double summation in (I80L i.e. the terms on the first 
line, and assume that off-diagonal terms average to zero due to the cosine function, then we find 


^ df / ~ . r,~ A 1 [mA/^ao (2m^-I-m3i) 


Vdt 


{-Px + 3Px) ) = E 


256Mpj7r 


1 - 


. 2 \ 1/2 
4mr, \ 


ruA 


(82) 


and by using the fact that pA = m^(Q;^) 2 / 2 d^, we can then express this as 


/ 1 ~ , o~ n\ V- 1 [/A(2m2 +TO^)]' 

\2/ dt^ V " E 128M4j^ 

i.e. we have found 


1 - 


,1 2 \ 1/2 
4mr \ 


rtiA 


m, 


Pa — ^ (83) 

A 


fA{^m\ + m\) 


a ->XX 


1287rMpjr7i^ 


1 _ Ek 
. ■ 


2 \ 1/2 


which is in agreement with (1361) . 


(84) 


The production of fermions 

In considering the case for fermions we follow closely the analyses of As many aspects of the 

calculation are similar to the bosonic case, we defer details to Appendix [BJ 
From the fermionic action in ([2]) we obtain the Dirac equation 

= 0. (85) 

Specialising to the case of a FLRW metric, and taking e° = 5° and e)j = Sl/a{t), we have 

7 °(x)= 7 °, y(x) = ^, Fo=0 and F, = ^7V- ( 86 ) 

a[t) Z 

If we introduce conformal time, and also define T = the Dirac equation can be written as 

{Ada + amp) T = 0, (87) 

where in this equation do = dr/. Correspondingly, the action and Hamiltonian can be written as 

5^ = — / dpd^x'i' {Ada + amp) 'k and = — f d^x^Ad^^ “ * / d^xAdrj'i’. ( 88 ) 
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The space of solutions is endowed with a conserved scalar product, and in the FLRW case it reduces to 

(4'i, 4'2) = y d^a:4'l4'2. (89) 

Given one solution to the Dirac equation, Ur(k,x), one can show that the charge conjugate Vr{k,x) = 
CU^ (k,x) = 7^17*(fc, a:), where C = 7^/3, is also a solution. Note that the subscript r labels the spin. We 
can then construct a basis of the solution space out of Ur(k, x) and Vr(k, x), and further require that the basis 
be orthonormal with respect to the above scalar product. As such, a general solution can be decomposed as 


4f(a:) = ^ f dPk(ar{k)Ur{k,x) + bl{k)Vr{k,x)^ 


(90) 


where ar{k) and 6)(fc) now correspond to annihilation and creation operators satisfying the anti-commutation 
relations 

^ar{k),al{^'^ = 6rsS''^'>{k - ^ and ^br{k),bl{^'^ = Srs6^^\k - (91) 

with all other commutators vanishing. We next decompose the solutions Ur{k,x) as 


Ur{k, x) = 


1 f UA{k,r])hr{k) \ a 


(92) 


(27r)3/2 \ rui3{k,r])hr{k) 
where k = k/k and hr{k) are the eigenvectors of the helicity operator 

k ■ ahrik) = rhr{k), r = ±1, (93) 

which are chosen to satisfy h\.(k)hs{k) = Srs- For the choice of hr(k) made in Appendix iBl we then find 

Imposing the orthonormality conditions dictates that 

|u>t(fc,?7)p-h |MB(fc,?7)P = 1, (95) 

and the Dirac equation now takes the form 


which can be decoupled to 


„ ( UA{k,r]) \ _ f am^ k \ f UA{k,f]) 
'^\UB{k,r])J k -am^ J \uB{k,r]) 


v) = - ± iiam^Y] UA,B{k, 77 ), 


(96) 


(97) 


These two equations are now of the same form as Eq. o for the boson mode functions. As such, the 
procedure from here onwards is very similar to the bosonic case. In analogy with the the bosonic case, we 
expand UA{k, rj) and ujsik, rj) in terms of positive and negative frequency functions as 


Wk am^ i j Wkdrj' 

2wk 


UA{k,v)= A, (77 )f - B, (77)y 

nB{k,v) = A(77)y^^^^e-*/“^-^’’' +Sfe(77)y^^^^e*/-^'^’'', 


(98) 

(99) 


where wf. = k'^ + a^mYp. With this decomposition we find that Ak{'rj) and Bk{ri) must satisfy the normalisation 
condition \Ak{'rii)\^ + \Bk{ji)\^ = I and the evolution equations 




( 100 ) 


( 101 ) 
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Assuming that at some time in the past Bk{r]Q) = 0, the above mode functions then coincide with the 
flat-space mode functions and the Hamiltonian is diagonal. As Bk{r]) evolve away from zero, however, the 
Hamiltonian is no longer diagonal, instead taking the form (49l - l5lj| 


H^ = ^ / (fk Wk (|Afc(?7)P - |Hfc(?7)p) (^al{k)arik) - br{k)bl{k)^ (102) 

-2Ak{r])Bk{r])wke~'‘'*’-f‘br{-k)ar{k) - 2A*k{r])Bl{r])wke^‘^-^a'^{k)bl{-k) 

In order to diagonalise the Hamiltonian one can make a Bogoliubov transformation, defining 

ar{k,r]) = Ak{ri)ar{k) - Bl{r])A^-i‘bl{-k), bl{k,r]) = Bk{'n)e~'’'*’kar{-k) + A*k{r])bl{k). (103) 


One then finds that the number operator associated with the new basis is given as {0\al{k,ri)ar{k,ri)\0) = 
{0\bl(k,r])br{k,r])\0) = |Hfe(? 7 )p. We must therefore determine Bkir]) if we wish to determine the number of 
particles created. Looking at (IIOII) we see that the form of the equation we need to solve for Bk{rf) is almost 
identical to that for Pkiv) that we solved in the case of the bosonic fleld0 As such, we defer details of the 
calculation to Appendix [B1 stating only the main results here. 

First, on using the stationary phase approximation we And that Bk{r]) is given as 


SfcW = +^SkT^/^] - Vo), 

A 


B 


A 

k 


^^piiwliVk) \ 


2T^a{Vk) 


(104) 

(105) 


where tpkiiv) is still as defined in (1571) but with = k"^ + a^rn^. As in the bosonic case, is the time at 
which tpk.'iiv) = 0 is satisfied, i.e. the time at which the phase is stationary, and is the sign of tpk.iiVk)- 
Then, as with the scalar case, we wish to determine the quantity {df /di)T^/{2f), which corresponds to the 
right-hand side of the continuity equation for the energy-momentum associated with ip in the Einstein frame, 
i.e. corresponds to the particle production term. Taking (0|T(’^)'',y|0) = diag(—pp, pp, pp), we find that 
to first order in Bkir]) 


= -pp + 3pp 


4 

(27r)3d^ 




(106) 


Proceeding in exactly the same way as for the bosonic field in the previous subsection, we arrive at 

1 df 


2fdi 


ifA)^rnlmA 

(-pp+3pp) 1 = 2 ^Pa -- - II- ' 


327rM4j 


from which we deduce 


Ua)^ 






327rM4j 


1 - 


4m: 


2 \ 3/2 


m2; 

A 


(107) 


(108) 


which is in agreement with (1361) . 


Energy-momentum tensors and their (non-)conservation 

In the analysis of the preceding two subsections we considered the continuity equation for the matter 
energy-momentum tensor in the Einstein frame. Our reason for doing so was the transparent interpretation: 


Indeed, if we had considered a conformally coupled field instead of a minimally coupled one, the evolution equation for the 
Bogoliubov coefficients would be the same up to a factor of fc/(ma) 14(1 . 
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the non-conservation of the matter energy-momentum tensor in the Einstein frame is a result of the explicit 
interaction terms that give rise to particle production. To close this section we consider the continuity 
equations for the other energy-momentum tensors. 

First let us consider the energy-momentum tensor of the oscillating fields in the Einstein frame. Combining 
the results of the previous two subsections we have 

= (109) 

,4 


As the total energy-momentum tensor must be conserved, this implies that aPa, where 

= Y.A^'^^si-PA^PA^PA^pA), with PA as given in (021) and 



( 110 ) 


Averaging over several oscillations we have {Pa) = 0, so that we obtain the expected continuity equation 
(HU). As previously mentioned, the interpretation of this standard result is intuitive - the energy density 
of the oscillating fields decays both as a result of the Hubble expansion and the decay into matter. Under 
the instant decay approximation, we assume that reheating ends once the decay rate “catches up” with the 
Hubble expansion, i.e. when = 3H. This then allows us to determine the reheating temperature in terms 
of Ta, and thus put constraints on model parameters such as Ja and niA- 

Next we consider the matter energy-momentum tensor in the Jordan frame, which we know to be covariantly 
conserved - recall (HU). At first glance this property would seem to be at odds with the fact that we have 
particle production. However, in the Jordan frame the particle production is interpreted as being due to the 
oscillatory nature of the Hubble rate, and the term on the right-hand side of dZSl), for example, becomes part 
of the Hubble dilution term on the left-hand side of the continuity equation in the Jordan frame. Namely, 
using py_ = , p-^ = and the relations given in (E21), dZSl) can be re-written as the standard continuity 

equation 


^^(«Vx) + H{-p^ + 3pp,) = 0. 


(Ill) 


Note that whilst we have considered the bosonic field as an example, the same is also true for any matter 
field. 

Finally we consider the energy-momentum tensor for the oscillating fields in the Jordan frame. In the 
Jordan frame there is some ambiguity as to how we might like to define the energy-momentum tensor of the 
inflaton fields, and the relation between the Jordan and Einstein frame inflaton energy-momentum tensors 
is not just a simple factor of as it is for the matter energy-momentum tensors. As commented in 
Appendix [21 Einstein’s equations in the Jordan frame can be recast into the standard form if we define the 
effective energy-momentum tensor given in (IA10|) . We then choose to define such that 


rp{eS) 




-^Pl rp(m) 


( 112 ) 


i.e. we have 




As such, we see that despite the fact that is covariantly conserved, is not, with 

A/r2 


In a FLRW background we explicitly have 

' 1 . 


P<p 


eff _ -^Pl '' 


/ 


+ V-3Hf 


eff _ -^Pl '' 


M, 


/ - nr'rdf 3 f df 


PI 


Kbr^’^-V + f + 2Hf 




dt 4/ \dt ^ 

r.rrdf 5 / rf/^ ^ 


(113) 

(114) 

(115) 

(116) 
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and the continuity equation 


pf + iH{pf+pf) = (117) 

where pf = pf = p^ = -f p^ = and The physical 

interpretation of this last equation is less clear than that of in the Einstein frame. However, one can 
assume that reheating completes when p^'^'> k. 3(p^® -\-pf)/2, where we have used the fact that f /f ^ ~277, 
as can be seen from (1221) . 

In the above analysis we have derived conditions for instant reheating in both the Jordan and Einstein 
frames. However, we note that imposing instant reheating in this class of models gives rise to issues regarding 
the discontin uity of H or H, which results from the assumption that / —>■ Mpj instantaneously at the time 
of reheating [l3|. To avoid this issue one must therefore solve the continuity equations dynamically. 


IV. SUMMARY AND CONCLUSIONS 


The high-precision nature of current CMB data dictates that reheating dynamics must be taken into 
account when trying to constrain different models of inflation. Given the recent interest in inflation models 
containing a non-minimal coupling to gravity and potentially multiple scalar fields, in this paper we have 
revisited the process of gravitational reheating that is inherent to this class of model. Our formulation 
allows for multiple, non-minimally coupled inflaton fields endowed with a non-flat field-space metric, and it 
is assumed that these Helds are not directly coupled to matter. 

At the level of the background dynamics, we saw that the oscillation of the inflaton fields about their vacuum 
expectation values gives rise to matter-dominated-like evolution of the Hubble rate in the Einstein frame, 
as in elementary reheating scenarios. In the Jordan frame, however, this matter-dominated-like evolution 
is modulated by an oscillatory component, and it is this oscillatory part that gives rise to the gravitational 
particle production of minimally-coupled matter, i.e. gravitational reheating. When interpreted in the 
Einstein frame the gravitational reheating does not result from the oscillatory nature of the Hubble rate, but 
instead from the explicit interaction terms between the inflaton sector and ordinary matter that are induced 
by the conformal transformation. 

In order to calculate the rate of particle production we used the method of QFT in a classical background, 
which requires the calculation of Bogoliubov coefficients. Although this was not entirely necessary for the 
perturbative reheating regime considered, the advantage is that much of the discussion will also carry over to 
the resonant preheating regime, where the perturbative flat-space QFT calculations are no longer applicable. 
Taking appropriate limits, we were able to confirm agreement between the Bogoliubov and perturbative 
QFT approaches, including kinematic suppression factors. Despite the difference in interpretation between 
the Jordan and Einstein frames, we saw that the calculation of the Bogoliubov coefficients associated with 
particle production was independent of the frame in which we started. This resulted from the fact that the 
canonically normalised quantum fields one naturally defines in the two frames are identical. 

To finish, let us mention one possible extension of the framework developed here. In analysing the dynamics 
of the oscillating inflaton fields at the end of inflation we made use of the mass eigen-basis of the Einstein 
frame potential. We implicitly made the assumption that all of the fields begin oscillating about their 
vacuum expectation values at approximately the same time, with mA ^ niB ^ H for all A and B. Such 
an approximation, however, may not be valid. Generally we might expect there to be a wide range of field 
masses, and that different fields therefore begin to oscillate and decay at different times. In the case that 
heavier fields are present, which start to oscillate and decay much earlier, it is perhaps reasonable to assume 
that the resulting decay products are diluted by inflation - which continues to be driven by the lighter fields 
- and are therefore negligible. However, in the case that lighter spectator fields are present, which do not 
oscillate and decay until much later, we might expect them to play a significant role. In general we would 
expect the field space metric hab, potential V{(j)) and non-minimal coupling function f{4>) to all depend on 
these spectator fields. Consequently, quantities such as S'ablvev and Vahivev, which we took to be constants in 
the analysis of Sec. IHI A[ would all become functions of the spectator fields. Ultimately, this would then lead 
to a spectator-field dependence of and through their dependence on niA and /a, which 

would in turn give rise to a modulated reheating scenario. We leave further consideration of this scenario to 
future work. 
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Appendix A: Energy-momentum tensors and their (non-)conservation 

In this appendix we review in more detail the relation between energy-momentum tensors defined in the 
Jordan and Einstein frames, including determining whether or not they are covariantly conserved. We follow 
closely [52l - [^ . simply generalising to the multi-field case. 


Covariant conservation of the matter energy-momentum tensor in the Jordan frame 


Let us start by showing that the energy-momentum tensor for matter in the Jordan frame is covariantly 
conserved, despite the presence of the non-minimal coupling. In Sec. IIIBi we derived the Einstein equations 
in the Jordan frame as 


G 




I 

7 





where, using the definition 




we have 


7 ^( 0 ) 







We also gave the equations of motion for the fields as 

- K + faR = 0, 


(Al) 


(A2) 


(A3) 


(A4) 


where F(jb|c = and F))^ is the Ghristoffel connection associated with the field-space metric hab- Taking 

the covariant divergence of (EH), and using the Bianchi identity, we find 


= 0 = \ + DV,/ - v,n/ 


(A5) 


Similarly, taking the covariant derivative of (EH and using the equations of motion (EH, we also find 




Substituting this result into (IA5I) . and recalling Gpu = R^i, — ^g^i^R, we arrive at 

v^Gp, = 0 = + J + DV,/ - v,n/ 

However, from the definition of the Riemann tensor we have 

□V,/ - v,n/ = 


which clearly then leaves us with 




= 0 = 


(A6) 


(A7) 


(AS) 


(A9) 
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i.e. we have recovered the fact that the matter energy-momentum tensor in the Jordan frame is covariantly 
conserved. 

One could also consider an effective energy-momentum tensor defined by jni as 


7^(eff) 

flL' 


Mi 

/ 



gij.uOf , 


(AlO) 


such that Einstein’s equations take the standard form This effective energy-momentum 

tensor is of course covariantly conserved as a result of the Bianchi identity. 


Energy-momentum tensors in the Einstein frame 


Turning to the Einstein frame, if we assume that only depends on and not its derivatives, then 
we can write 

_2 J _2 J 

SgP"^ dgP''' 

Under the conformal transformation we have g^i, = g^'^ = g*^'^ and \J—g = which we can 

substitute into (lAllIl to find 

= (A12) 

Recall that = f{cf))/Mp^. Taking the covariant divergence of this (with respect to the Jordan frame 
metric) we have 


y/xjd™) = . (A13) 

We now need to use the relation between covariant derivatives as defined with respect to gf^, and The 
relation between the Christoffel symbols is given as 

i , (AM) 

and on substituting this result into ()A13ll we hnd 

^ Q ^ -p ^ (ai5) 

where tI™! is the trace of the matter energy-momentum tensor. Thus, we see that even if is covariantly 
conserved, in general is not. It will, however, be conserved if tI™) = 0, which is the case for radiation-like 
matter. 

Given that in the Einstein frame we have 

V'^G^, = 0 = (rW + , (A16) 

the non-conservation of implies a non-conservation of . Let us try to determine this explicitly. 

When we try to calculate the equations of motion for the scalar fields in the Einstein frame, we need to 
correctly take into account the dependence of y/—gL^"^'> on (j)^ that results from the conformal transformation. 
However, as the only dependence of yf—gC^'^^ on </>“ comes from the conformal transformation 
if we once again assume that y/—gC^'^'> only depends on g^i, and not its derivatives, then we can simply use 
the rules of partial differentiation to obtain 


J ^ J dgP'' ^ Sjy^gC^ ^ 

SgP^ 5gp‘' fl ^ SgP’' 

As such, the equations of motion for the scalar fields become 


(A17) 


- SabDcf» + Ua + = 0 


n 


(A18) 
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where is the Christoffel connection associated with Sab- As given in the main 

text, the energy momentum tensor T^t'^ takes the form 

= SabV^r^^cj)^ - Qfaa V,,/ + . (A19) 

Taking the covariant divergence of this and making use of the equations of motion (IA18I) . we find 

(A20) 

This is entirely consistent with what we expected from (jA16l) and (IA15I) . 


Appendix B: Details regarding fermion particle production 

In this appendix we give additional details regar ding the fermion particle production calculation. As stated 
in the main text, we follow closely the analyses of 


Conventions 


First let us clarify our conventions, 
representation 


We use the (—h ++) sign convention and the Dirac gamma matrix 



(Bl) 


where cr® are the Pauli matrices. 

Following , for the eigenvectors of the helicity operator we use the following spherical-coordinates-based 
representation 


h+i{k) 



h-i{k) 



(B2) 


where we have chosen the normalisation such that h].(k)hs(k) = 5rs is satisfied. With this choice, it is then 
possible to show that 

-ia^h*{k) = -rC*-’^h_r{k), (B3) 

which on using hr{—k) = h-r{k) allows us to find 


Vr{k)=J^U:{k,x) 


f -u*^{k,ri)hr{-k) 

(27r)3/2 ru\{k,'q)hr[-k) 


-ik-x 


(B4) 


Calculation of production rate 

Here we give more details regarding the calculation of the production rate for fermions. As with the 
scalar case, in order to determine the particle production rate we need to evaluate the right-hand side of the 
continuity equation for the fermion energy-momentum tensor in the Einstein frame. As such, we need to 
first determine the vacuum expectation values of the components of the energy- momentum tensor, which we 
denote (0|f= diag( —Pi/j, PtP, Pip, Pip)- In terms of the quantities pip and pip the continuity equation 
takes the form 


Pip) + H\~Pip + “ipip) = A ^Pip)- 


2/dr 


(B5) 
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The explicit form of for fermions is given in (I20L and on writing in terms of the canonically normalised 

field d'(x) given in (1901) and taking the vacuum expectation value we find 


Pip 

Pip T ^Pip 


(27r)3 


d^kwk \}3kir])\ 



idm^ 

d'^fl 


(Pk am^ /1 
(27r)3 _Wk^ V2 



Wk ^ / 


(B6) 

(B7) 


If we are assuming that initially no fermion particles are present, then it is appropriate to consider the 
perturbative regime where Bk{p) -C 1 and Ak{p) ~ 1. As such, to leading order in Bk{p) and we find 


^ df / ~ o ~ N ^ df 4(47r) 

—+ 3p^) 


'dkk^^^?fi 

Wk 




—2i f Wkdrj' 


(B8) 


It is now clear we must solve for Bk{p), i.e. solve Eq. (IIOII) . In order to do so, we note that to leading 
order in cP we have 




.A' 


( V ~ f fj 

[amA -am^ „ + a Um^. 


2M2j 


(B9) 

(BIO) 


Note that, as with the bosonic case, taking expressions to leading order in ensures that we are only 
considering the tri-linear interaction terms and the perturbative regime appropriate for comparison with the 
perturbative QFT calculations of Sec. IIIIBI Taking Ak{p) -A I and integrating (IIOII) we have 


Bk{v) 



~2£V2 ~fA<A 

an —a 


Ai 


2M2j 


c-'^Ad^Wkdri" 


(Bll) 


The first term in the brackets is slowly varying, so that its contribution to Bk{p) averages to zero. The second 
term, however, is highly oscillatory, thus giving a non-zero contribution to Bk{p) that can be calculated using 
the stationary phase approximation. Explicitly, we have 


Bkiv) 


E 

A 


r ^ , kpAf^rnAaom^ 
Im 8A/2jm2 


JiriAipkAv') _ JmAiptAn 




(B12) 


where '4’ki ' 4 ’k 2 given in (Ell) and (1681) but with wf. = P + a^rn^. The phase no stationary 

points for physical values of a, and therefore the second term in the above expression gives no contribution 
to Bk{p)- The first term involving the phase however, does have a stationary point at a = 2wk/mA- To 
leading order in a^, this stationary phase condition gives 


k niA f 


(B13) 


where rj^ denotes the time at which the condition is satisfied for a given k and niA- As in the scalar case, 
this coincides with our expectation from kinematics. In making the stationary phase approximation one also 
needs the second derivative of the phase at the time 77 ^, which in the fermionic case is given as 



a^{Pk)H{pt) 



+E 

B 


ml 


(B14) 


The final solution for Bk(rj) is given in (110411 . 

We next turn to evaluating (IBSp . On expanding df /dt we obtain 


1 df _ 

2fdi^ 


■ ^Pip) — 


4(47r) 

(27r)^Q 


oIVl-p^l'W L \ 


(B15) 
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We see that this is a highly oscillatory function, and as with the scalar case we consider taking an average 
over several oscillations as 


1 df . 
2fdi^ 


3P^) ) = ^ / ^ 


rv+T 


lr,-T 


4(47r) 

(27r)3d^ 




X 2iQ 




(B16) 


where T ^ 0(l/(dmA))- This averaged quantity will only be non-zero if rj coincides with a stationary point 
of the phases ipkiiv) or '4’k2iv)- Seeing as V'fe 2 (^) >^0 stationary points for physical values of d, we only 

obtain contributions from terms involving On making the stationary phase approximation we arrive 

at 


2fdi 


(-Pb + 3p ^)) =Y1 


4(47r) 


A,B 


(271)^0^(77) 


dkS{r] - rii)wkiVk)k^ 


(B17) 


X 25R 


)+2(sfc-sf )7r/4 


k 

X 0(p^-pf)0(pf -Po), 


where is as defined in (IlHhIl . On using the fact that = B^*B^, this can then be written as 


1 df , 
2fdi^ 


3p^) 


= E J ~ dk)Mrik)k'^\Bt? 

X 2 cos {mAf’tii'nk) - ^BipkAijlk) + (Sfe - sf )7r/4) 6 ( 77 ^ 


(B18) 


Po)- 


In the above expression we have arranged that ms > bia for B > A, meaning that 6 ( 77 ^ ~ Vk) — ^ only for 
B > A. Next we re-write the delta function in 77 as a delta function in k. The relation is as given in (ISTl) . 
but now with tp^iiVk) given in (IB14I1 and ija to leading order in given as 



(B19) 


Assuming that the non-diagonal terms in the second line of (IB18|) average to zero, the diagonal terms give 
rise to (I107D . 
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